Optical conductivity near finite-wavelength quantum criticality 
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We study the optical conductivity a (SI) of an electron system near a quantum critical point with 
finite-wavelength ordering, a (SI) vanishes in clean Galilean-invariant systems, unless electrons are 
coupled to dynamical collective modes, which dissipate the current. This coupling introduces a non- 
universal energy scale. Depending on the parameters of each specific system, a variety of responses 
arise near criticality: scaling peaks at a temperature- and doping-dependent frequency, peaks at a 
fixed frequency, or no peaks to be associated with criticality. Therefore the lack of scaling in the 
far-infrared conductivity in cuprates does not necessarily call for new concepts of quantum criticality. 

PACS numbers: 74.25.Gz, 71.45.Lr, 73.20.Mf, 73.43.Nq 
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Quantum criticality has become a central issue in solid- 
state physics and might account for the non-Fermi-liquid 
properties of cupratesi^i^, heavy fermions^ and other 
strongly correlated electron systems. A relevant goal is 
to determine the role of critical fluctuations in various 
anomalous properties. In this context, controversial re- 
sults on the scaling behavior of optical conductivity a(Sl) 
have been found^^ i 9 ' 10 ! 11 ] 12 . The issues of scaling viola- 
tions and the possible need of "concepts beyond the stan- 
dard model of quantum criticality"— have been discussed 
on the basis both of experimental and theoretical ar- 
guments. In this Rapid Communication we consider the 
optical conductivity of a system near a quantum critical 
point (QCP). We show that difficulties in reconciling op- 
tical conductivity experiments with the standard quan- 
tum criticality framework does not necessarily call for 
new concepts of quantum criticality, but may naturally 
arise from the specific nature of this response function: 
a(Q) measures absorption and requires a mechanism for 
current dissipation (impurities, umklapp processes, dy- 
namical phonons all play this role in real systems). This 
mechanism, in turn, involves additional energy scales 
possibly leading to violation of scaling in some param- 
eter range. 

Our analysis is carried out for a neutral order pa- 
rameter with critical collective modes (CM's) at finite 
wave vectors q c . This is not only pertinent to the 
cuprates, which are the main object of this work and 
where a relevant role of spin-orderingi and/or charge- 
ordering 2 fluctuations has been proposed, but may apply 
to some heavy fermions, dichalcogenides, and other sys- 
tems where spatial (usually spin or charge) order occurs 
at low temperature. The neutral order parameter fluc- 
tuations couple in a non-trivial dynamical way with the 
external electromagnetic field. 

We address two issues: (i) to what extent critical CM's, 
usually characterized by a strong dependence on tem- 
perature and other control parameters (like pressure or 
chemical doping x), produce specific signatures in a(fl) 
and (ii) how dissipation mechanisms introduce additional 
scales, which may mask the quantum critical features of 
optical spectra in some parameter range. 

In the first part of this Rapid Communication we con- 



sider the textbook case of Galilean-invariant electrons, 
where current dissipation coincides with momentum dis- 
sipation. We will later see that the Ward identities 
enforcing current conservation are in fact more general 
when implemented for the leading contribution arising 
from scattering due to CM's at finite q c . In this case our 
results hold true for electrons in a lattice whenever the 
umklapp processes for critical scattering are inactive. 

- The model — We start from a model effective action 
at temperature T 

s = rj^w^c^ + rjxo" 1 !?)^-, 

k.a q 
q k,(7,a' 

where k = (k, ee), q = (q, u> n ), e£ (ui n ) are fermionic 
(bosonic) Matsubara frequencies, Go = («Q — Ck) _1 is the 
bare quasiparticle (QP) propagator, £t is the dispersion 
of the QP's created by the c fields, and 7 is the (theory- 
dependent) vertex coupling the QP's to collective charge 
or spin excitations, represented by the bosonic fields <fi. 
Xo is the bare CM propagator, i.e., the charge or spin sus- 
ceptibility in the absence of QP polarization dressing. At 
the bare level the <j> fields mediate an effective electron- 
electron interaction V(q) = 7 2 Xo(?)- If this interaction is 
purely static, it cannot lead to electron momentum dis- 
sipation in a Galilean-invariant system and, as we show 
below, o~(Q) vanishes at any finite frequency, provided 
it is properly calculated within a conserving perturba- 
tive scheme. If instead this bare interaction has its own 
dynamics arising, e.g., from phonons^^, electron mo- 
mentum is dissipated and the optical response is finite. 
For the sake of concreteness we consider the charge case, 
where the <fi fields and the 7 vertex are scalar. To make 
contact with the case of dichalcogenides and cuprates, 
we consider a model in two dimensions. We perform a 
conserving calculation starting from the simplest bubble- 
type Baym-Kadanoff (BK) functional 14 , which involves 
only bare QP propagators and the bare interaction V(q) 
[see Fig. QJa)]. Once current vertices are inserted, to cal- 
culate <r(f2) by means of the Kubo formula, the bubble 
diagrams can be resummed by introducing a random- 
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FIG. 1: (a) Typical diagram of our BK generating functional. 
The solid and dashed lines represent the QP propagator and 
the bare interaction V(q), respectively, (b) Bubble resumma- 
tion to obtain the RPA-dressed CM propagator (wavy line), 
(c) Diagrams for the current-current response, obtained from 
current-vertex insertions (solid dots) in the diagrams of the 
BK functional. AL, V, and S stand for Aslamazov-Larkin- 
like, vertex, and self-energy diagrams, respectively. 



phase-approximation- (RPA) dressed CM propagator [see 
Fig. Hfb)]. Then, the full set of conserving diagrams re- 
ported in Fig. QJc) is obtained. These diagrams give the 
CM corrections to the bare current-current response. In a 
Galilean-invariant system, er(il) maintains a <5-like Drude 
term, while the critical CM's may give absorption at fi- 
nite frequencies. 

The RPA effective interaction of Fig. [TJb), T = (V^ 1 - 
LTr 1 [here 11(g) = -T£ fc G (fc + q)G (k) is the QP 
polarization bubble] , near the QCP has the generic form 



T(q) = -TlZ 1 (m + e q 



(i) 



for small frequencies and q ±q c - Here, LT W = 
[II(q c ,cJ„) - II(q c ,0)]/ \u n \ | Wn =o, and e q « z/|q =F q c | 2 . 
to, v, and Vl are model-dependent coefficients which re- 
sult from the second-order expansion of V~ Y — II around 
q c and u n = 0. Except for the uj 2 term, T(q) has 
the general hydrodynamic form of a propagator for dif- 
fusive CM's (damped by QP's) near a Gaussian QCP. 
Here to, proportional to the square of the inverse cor- 
relation length, is the CM mass measuring the distance 
from criticality. If the frequency dependence of T only 
arises from the QP bubble n, i.e., if the bare interac- 
tion V is purely static, the optical response is zero [the 
diagrams of Fig. 1(c) cancel each other], as we show 
below. On the other hand, if the bare interaction V 
has its own dynamics, a finite <r(f2) is obtained. We 
analyze these two situations, and for concreteness we 
fix the values of the parameters, adopting as an ex- 
ample the model of Ref. [TH There, a bare interaction 
V(q) = Vb(q) — \ui 2 (ui 2 + oj 2 ) -1 was considered, arising 
from both a static short- and long-range Coulomb repul- 
sion Vq, and from the coupling to a dispersionless phonon 
of frequency Co. For moderate electron-phonon coupling 
A ^ £ f, where £p is the typical QP Fermi energy (in 
cuprates, e.g., Ef ~ 0.3 eV), this interaction can lead to 
a charge-ordering instability, at a wave vector q c . Within 
this model we find to = n" x {[A - Vb(q c )] _1 + H(q c ,0)}, 




FIG. 2: Sketch of a two-dimensional Fermi surface with four 
HS's connected by a generic incommensurate critical wave 
vector q c = (q c ,0). Fermi velocities at the HS's and their 
components are also shown. 



and fi = A" 1 ^ - V (q c )] 2 n w . In Refs. and [H the 
instability, signaled by a vanishing to, occurred for A ~ 
Vb ~ |TI(q c , 0) | 1 ~ Ep. Then, we estimate v ~ Epjk 2 F 
{hp is the Fermi momentum), and f2 ~ uj(lu/ef). 

- Current- current response function — Near crit- 
icality, the diagrams of Fig. QJc), with incoming 
zero momentum and finite frequency f2;, are domi- 
nated by the poles of the CM propagators. Then, 
we write the vertex-self-energy (VS) and Aslamazov- 
Larkin (AL) contributions to the current-current re- 
sponse function x"° as T S q , w „ Vyg(o;„, Oz)r(q, w n ) 
and |TE q , Wn [V5 L K,O i )] 2 r(q,'^)r(q,n / +a; n ), where 
a = x,y and we have exploited the relation r(q, u> n ) — 
r(— q, u> n ). The vertices Vyg and V% L come from the 
integration on the QP loops. To perform an analytic cal- 
culation, we adopt the standard procedure of linearizing 
the QP dispersion around the points of the Fermi surface 
connected by q c [hot spots (HS's)Ji. Then, we find (see 
also Ref. GJ) 



V£S(w n> n,) 



x (\Qi+uj n \ -2|w w | + |fi,-w„|), (2) 
VgtK.n,) = -ieD^^nf^lfij+Wnl-lwnl), (3) 



where e is the electron charge, u a 



— U HS1 



J HS2 



, and 



vfj S is the a component of the Fermi velocity at the 
HS's (see Fig. [5]). The vertices vanish identically in the 
direction perpendicular to q c , i.e., for a = y, since in 
Fig. [5]we took q c along the x axis and v v HS1 = v v HS21 i.e., 
U v = 0. Thus Xtfitil) = for fl t ^ 0, regardless of the 
retarded or static character of the bare interaction V. 
To calculate Xjj> we exploit the identity 



r(q,w n )r(q, Qi+uj n 



nj 1 [r(q,a+^ 1 )-r(q,q; ?t )] 

\n l + u n \-\u n \ + ^(n l + 2u n y 

(4) 

which is immediately derived from Eq. |T]) , and allows us 
to write the AL contribution to Xjf as 



|0;+o; n | - \w n \ + 2t(ni + 2u n ) 



r(q,w„), 



where a single effective interaction appears, at the ex- 
pense of a more complicated prefactor. 

It is instructive to consider first the case Q — ► oo, 
(i.e., u> — > oo). In this case the QP's are only coupled 
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among themselves and with a static phonon held, and 
the frequency dependence in Eq. |T]) only comes from 
the QP bubble. As a consequence, <j(Sl) must vanish for 
any finite frequency, while only the 5-like Drude term re- 
mains. For SI — ► oo and by means of Eq. ([3]), Eq. ((4]) 
can be cast in the form of a Ward identity eu"[T(q, Sli + 
w n ) -r(q,w„)] = i^liV% L {uj n ,Vli)T{ci, Sl t +w„)r(q,w„), 
which plays the role of the ordinary Ward identity ex- 
ploited, e.g., in Ref. [171, to prove the cancellation of 
the paraconductivity, in the case when the electron mo- 
mentum commutes with the Hamiltonian. In our case, 
by means of Eqs. @ and © we can eliminate the 
VS vertex in favor of the AL vertex, Vy x s (ui n ,Sli) — 

-^« B fir 1 [VALK,ni) + ^li(-wn,n I )]. Then, x j; = 
fu^TE^jv^K,^) - v^(-^,^)]r(q,^) 

vanishes identically, since the quantity between square 
brackets is an odd function of uo n . 

- Momentum dissipation: finite SI — To obtain an 
optical response it is therefore necessary to introduce a 
dissipation mechanism. In this Rapid Communication 
we implement this dissipation via the w 2 term in Eq. |T]), 
which in the example we are considering arises from the 
phonon dynamics. In this case, after analytic continua- 
tion of the external Matsubara frequency, we find 



A 



, +oc dz 
PI —In 

ITT 

— oo 



(iSl — Sl)Cl J_oo tir 

Usi + i(2z + n)(z - n) 
n-i(2z + fi) 

'g* -"><*-"> coth (*— 
Sl-i(2z-Q) \ 2T 



z 2 — n(m 
coth (—) 

si 



(5) 



where A = e 2 (vjj S ) 2 jvd is the dimensional prefactor, 
d is the interlayer distance which translates the two- 
dimensional response into the in-plane response of a lay- 
ered system, and A is an ultraviolet cutoff. Extracting 
the conductance quantum e 2 /h and the factor 1/d, we 
are left with a dimensional factor which we estimate as 

Kv x hs) 2 /v ~ £F/h. 

One can check that Im Xjj (SI) linearly vanishes with SI, 
giving a finite cr(f2) = lmXjj(Sl)/Sl, for SI — ► 0. However, 
in our calculation without disorder, a 5-like Drude term 
[Dq — nReXjj(0)]S(Sl) is still present, where D is the 
QP Drude weight in the absence of the CM contribution. 
Kramers- Kronig relations connecting Re Xjf and Im Xjj 
guarantee spectral weight conservation within our con- 
serving approach: the finite-frequency weight associated 
with a(Sl) is exactly subtracted from D Q by 7rRexJ?(0), 
i.e., / dSl {[Do - 7rRe X J/(0)] S(Sl) + a(Sl)} = D . 

In Eq. (|5|), the SI dependence is crucial and introduces 
a non-critical energy scale in the absorption, which may 
strongly alter the dependence of the optical spectra on 
the CM mass. The scale SI (which we estimated for the 
cuprates as a rather low energy scale ~ u) 2 /ef) deter- 
mines, at low T, the frequency below which the response 
tends to be vanishingly small, as it is for SI = oo. On the 
other hand, for SI > SI a finite absorption is found. We 
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FIG. 3: (Color online) Optical conductivity for a system of 
QP's coupled to CM's in the quantum-critical regime (x = 
xqcp), with m = 30T and ft = 30 cm -1 . All frequencies and 
temperatures are in cm -1 . Fhe interlayer distance was taken 
as d — 10 -9 m. 
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FIG. 4: (Color online) Same as Fig. [3] but with m = 3T and 
f2 = 330 cm -1 . 



discuss a system in the quantum critical regime where 
the mass m(T,XQCp) = m{T) = aT^ 

We first consider the case a = 30 and SI = 30 cm -1 , 
illustrated in Fig. [31 As long as m(T) > SI, the peak 
position is ruled by the critical energy scale m and the 
strong T dependence of the CM mass is visible as peaks 
in a (SI) at temperature-dependent frequencies: the peaks 
soften and sharpen upon decreasing T. Thermal effects 
tend to fill the spectra below the peaks, which are clearly 
visible only if the CM mass, i.e., the peak frequency, is 
substantially larger than T. We find indeed absorption 
peaks at frequencies SI ~ m(T) only if the non- universal 
coefficient a > 10, as in Fig. [3j On the other hand, 
in Fig. H we take m(T) = 3T < H « 330 cm" 1 and 
the peak position is fixed at SI.— In this case, thermal 
effects show up for T > 0.2S1 with absorption at low 
frequency. Then, the peaks are gradually embedded in 
this thermally-generated absorption. This shows that the 
absence of well-separated absorption peaks in the far in- 
frared optical spectra does not necessarily mean that crit- 
ical CM's contributing to a (SI) are absent. 

- Discussion — So far we considered clean Galilean- 
invariant systems. The presence of a lattice introduces 



4 



two effects: the current vertex is no longer proportional 
to the electron momentum and umklapp scattering pro- 
cesses may occur. At leading order in the critical scat- 
tering, the first effect is immaterial because the current 
vertex is evaluated at the HS's and the velocity only ap- 
pears in multiplicative prcfactors. Then, for Cl — > oo 
the Ward identities implementing the cancellation remain 
valid, while Eq. holds for finite Cl. The second effect 
may instead be significant. However, our results for the 
leading critical contributions remain valid in those sys- 
tems where q c cannot connect Fermi surfaces in different 
Brillouin zones. Our analysis also provides reliable re- 
sults for dirty systems in the frequency range above the 
typical impurity scattering rate 1/r. In assessing the re- 
sults discussed above, we strongly relied on a conserving 
scheme to calculate <r(Cl). Although the relevance of ver- 
tex corrections was previously proposed 20 , in this Rapid 
Communication we provide quantitative evidence that a 
conserving scheme captures the strong cancellations oc- 
curring in systems with current conservation and it keeps 
partial cancellations also when the current conservation 
law is not obeyed. Therefore, non-conserving approaches 
in models with strong momentum dependence 2 - should 
be handled with care^ 

Our main result is that (the leading critical contribu- 
tion to) optical absorption at finite frequency only ap- 
pears, in a clean Galilean-invariant system or in systems 
without critical umklapp scattering, when electrons are 
coupled to other dynamical degrees of freedom, which 
mediate a retarded interaction allowing for electron- 
current dissipation. This introduces a non-universal en- 
ergy scale Cl. If the CM mass to, controlling the dis- 
tance from criticality, is larger than Cl, the CM peak po- 
sition is ruled by to and follows its x and T dependence. 
Otherwise, the peak is pinned at a frequency Cl. This 
saturation, which occurs even when to vanishes with T 



(x = xqcp), is an intrinsic feature irrespective of other 
mechanisms (e.g., pinning), which would keep to finite. 
Therefore, it is crucial to recognize that the saturation 
of the peak in far infrared spectra cannot generically be 
used as evidence against criticality. Thermal filling re- 
duces the visibility of the CM peak to the regimes where 
the peak frequency is substantially larger than T. When- 
ever the above conditions are not fulfilled, the CM's pro- 
vide broad absorption. In real systems, such absorption 
would interplay with the Drude peak and could be con- 
fused with an additional impurity broadening. 

Experiments in cuprates provide a wealth of differ- 
ent behaviors, which may find counterparts in the var- 
ious regimes discussed above. In Bi 2 Sr 2 Cu06 (Ref. 0) 
and Nd2_ x Ce x Cu04_ y (Ref. H) peaks at temperature- 
dependent frequencies, clearly separated from the Drude 
peak, are observed, with approximate scaling behav- 
ior. Within our theory, this implies to > Cl. Peaks 
at temperature-independent frequencies (but still dis- 
tinct from the Drude peak) are found in La 2 Sr x Cu04 
(Ref. H), Lai,6-xNdo.4Sr x Cu04 (Ref. 03), and some 
YBa 2 Cu307_ y (Ref. Ill)) samples. This behavior might 
be due to a non-vanishing mass or to a saturation at Cl. 
No peaks at all are observed in Bi 2 Sr 2 Cu 2 0g (Ref.0) and 
YBa 2 Cu307_ y (Ref. 12). Of course, an explicit quantita- 
tive analysis of the optical conductivity at low frequencies 
requires a proper treatment of disorder and umklapp pro- 
cesses, which is presently in progress. Nevertheless our 
theory accounts for the variety of experimental behaviors 
and may reconcile standard quantum criticality with the 
lack of universal scaling of the conductivity. 
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lani and J. Lorenzana and financial support from the 
MIUR-PRIN 2005, prot. 2005022492 and from the 
Alexander von Humboldt foundation. 
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